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Abstract 

This article describes the REDUCE package ZEILBERG implemented by Gregor Stolting and 
the author. 

The REDUCE package ZEILBERG is a careful implementation of the GosperQ and Zeilberger 
algorithms for indefinite, and definite summation of hypergeometric terms, respectively. An 
expression is called a hypergeometric term (or closed form), if ak/cLk-i is a rational function 
with respect to k. Typical hypergeometric terms are ratios of products of powers, factorials, 
r function terms, binomial coefficients, and shifted factorials (Pochhammer symbols) that are 
integer-linear in their arguments. 

1 Gosper Algorithm 

The Gosper algorithm [|^ is a decision procedure, that decides by algebraic calculations 
whether a given hypergeometric term has a hypergeometric term antidifference gk, i. e. 
dk — 9k-i = Ofc, and returns Qk if the procedure is successful, in which case we call Gosper- 
summable. Otherwise no hypergeometric term antidifference exists. Therefore if the Gosper 
algorithm does not return a closed form solution, it has proved that no such solution exists, 
an information that may be quite useful and important. The Gosper algorithm is the discrete 
analogue of the Risch algorithm for integration in terms of elementary functions. 
Any antidifference is uniquely determined up to a constant, and is denoted by 

Finding given is called indefinite summation. The antidifference operator E is the inverse 
of the downward difference operator Va^ = at—ak-i. There is an analogous summation theory 
corresponding to the upward difference operator Aa^ = 0^+1 — a^. 
In case, an antidifference gk of is known, any sum 

n 

^ Ofc = gn — g-m-l 
k=m 

^Thc sum package contains also a partial implementation of the Gosper algorithm. 
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can be easily calculated by an evaluation of g at the boundary points like in the integration 
case. Note, however, that the sum 

sC) 

e. g. is not of this type as the summand (j^ depends on the upper boundary point n explicitly. 
This is an example of a definite sum that we consider in the next section. 
Our package supports the input of powers (a"k), factorials (factorial (k)), F function terms 
(gamma (a)), binomial coefficients (binomial (n,k)), shifted factorials (pochhammer(a,k) = 
a(a + 1) ■ ■ ■ (a + — 1) = T{a + k)/T{a)), and partially products (prod(f ,k,kl ,k2)). It 
takes care of the necessary simplifications, and therefore (in principle) provides you with the 
solution of the decision problem as long as the memory or time requirements are not too high 
for the computer used. 



2 Zeilberger Algorithm 

The (fast) Zeilberger algorithm deals with the definite summation of hypergeometric 

terms. Zeilberger's paradigm is to find (and return) a linear homogeneous recurrence equation 
with polynomial coefficients (called holonomic equation) for an infinite sum 

oo 

= J2 /(^'^) ) 

fc=— oo 

the summation to be understood over all integers k, if f{n, k) is a hypergeometric term with 
respect to both k and n. The existence of a holonomic recurrence equation for s{n) is then 
generally guaranteed. 

If one is lucky, and the resulting recurrence equation is of first order 

p{n) s{n — 1) + q{n) s{n) = {p,q polynomials) , 

s{n) turns out to be a hypergeometric term, and a closed form solution can be easily estab- 
lished using a suitable initial value, and is represented by a ratio of Pochhammer or F function 
terms if the polynomials p, and q can be factored. 

Zeilberger's algorithm does not guarantee to find the holonomic equation of lowest order, but 
often it does. 

If the resulting recurrence equation has order larger than one, this information can be used 
for identification purposes: Any other expression satisfying the same recurrence equation, and 
the same initial values, represents the same function. 

Note that a definite sum J2 /(^; k) is an infinite sum if f{n, k) = for k < mi and k > m2. 

k=m\ 

This is often the case, an example of which is the sum (|l]) considered above, for which the 
hypergeometric recurrence equation 2s{n—l)—s{n) = is generated by Zeilberger's algorithm, 
leading to the closed form solution s{n) = 2". 

Definite summation is trivial if the corresponding indefinite sum is Gosper-summable analo- 
gously to the fact that definite integration is trivial as soon as an elementary antiderivative 
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is known. If this is not the case, the situation is much more difficult, and it is therefore quite 
remarkable and non-obvious that Zeilberger's method is just a clever application of Gosper's 
algorithm. 

Our implementation is mainly based on [Q. Many more examples can be found in 0], [0, 
and [§] most of which are contained in the test ffie zeilberg.tst. 



3 REDUCE operator GOSPER 

The ZEILBERG package must be loaded by: 
1: load zeilberg; 

The gosper operator is an implementation of the Gosper algorithm. 

• gosper (f ,k) determines a closed form antidifference. If it does not return a closed form 
solution, then a closed form solution does not exist. 

• gosper (f ,k,m,n) determines 

n 
k=m 

using Gosper's algorithm. This is only successful if Gosper's algorithm applies. 
Example: 

2 : gosper ( (-1) " (k+1) * (4*k+l) *f actorial (2*k) / 

(factorial (k)*4"k*(2*k-l)*factorial(k+l)) ,k) ; 



k 

- ( - 1) *f actorial (2*k) 



k 

4 *f actorial (k + 1) *f actorial (k) 

This solves a problem given in SIAM Review (g]. Problem 94-2) where it was asked to 
determine the infinite sum 

o c c _^ (-l)^'+^(4fc + l)(2fc-l)!! 

^-i™^- ^--Z. 2^^2k-l){k+l)\ ' 

{{2k — 1)!! = 1 • 3 ■ ■ ■ {2k — 1) = -^r^). The above calculation shows that the summand is 

Gosper-summable, and the limit S* = 1 is easily established using Stirling's formula. 

The implementation solves further deep and difficult problems some examples of which are:0 



^Note that REDUCE Version 3.5 gives the output in terms of F functions. 
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3 : gosper ( sub (n=n+l , binomial (n , k) "2/binomial (2*n , n) ) - 
binomial (n,k) "2/binomial(2*n,n) ,k) ; 



2 2 
((binomial(n + l,k) *binomial(2*n,n) - binomial (2* (n + 1) ,n + l)*binomial(n,k) ) 

2 

*(2*k - 3*n - l)*(k - n - 1) )/( 

2 

(2*(2*(n + 1) - k)*(2*n + l)*k - (3*n + l)*(n + 1) ) 
*binomial(2*(n + 1) ,n + l)*binomial(2*n,n) ) 
4: gosper (binomial (k, n) ,k) ; 
(k + l)*binomial(k,n) 



n + 1 

5 : gosper ( (-25+15*k+18*k~2-2*k"3-k~4) / 

(-23+479*k+613*k~2+137*k-3+53*k-4+5*k-5+k-6) ,k) ; 

2 

- (2*k - 15*k + 8)*k 



3 2 
23* (k + 4*k + 27*k + 23) 



The Gosper algorithm is not able to give antidifferences depending on the harmonic numbers 

k ^ 

Hk-=Y.-^ 

j=iJ 

e- g- J2kHk = {k + l){Hk+i — 1), but, is able to give a proof, instead, for the fact that Hk 
does not possess a closed form evaluation: 

6: gosper (1/k, k) ; 

***** Gosper algorithm: no closed form solution exists 

The following code gives the solution to a summation problem proposed in Gosper's original 
paper [Q. Let 

k k 

/fc = n (« + + cj') and gk = H + + cf) • 

i=i i=i 

Then a closed form solution for 



fk-l 

k 



9k 

is found by the definitions 



4 



7: operator ff,gg$ 



8: let {ff(~k+~m) => f f (k+m-1) * (c* (k+m) ~2+b* (k+ni)+a) when (fixp(m) and ni>0) , 
ff(~k+~m) => ff (k+m+l)/(c*(k+m+l)~2+b*(k+m+l)+a) when (fixp(ni) and ni<0)>$ 

9: let {gg(~k+~m) => gg(k+ni-l)*(c*(k+m)"2+b*(k+ni)+e) when (fixp(m) and m>0) , 
gg(~k+~ni) => gg(k+m+l)/(c*(k+m+l)~2+b*(k+m+l)+e) when (fixp(m) and m<0)>$ 

and the calculation 

10: gosperCff (k-l)/gg(k) ,k) ; 

ff (k) 

(a - e)*gg(k) 

11: clear ff,gg$ 

Similarly closed form solutions of J2k ^^^f^ positive integers m can be obtained, as well as 
of Efc — for 

9k 

k k 

fk = l[('' + bj + cf + df) and gk = l[((' + bj + cf + df) 

i=i i=i 

and for analogous expressions of higher degree polynomials. 

4 REDUCE operator SUMRECURSION 

The sumrecursion operator is an implementation of the (fast) Zeilberger algorithm. 

oo 

• sumrecursion (f ,k,n) determines a holonomic recurrence equation for J2 f{^ik) 

k=—oo 

with respect to n. 

• sumrecursion(f ,k,n, j) (j G IN) searches only for a holonomic recurrence equation of 
order j. 

A simple example deals with Equation (P0 
12: suinrecursion(binomial(n,k) ,k,n) ; 

2*suin(n - 1) - sum(n) 

The whole hyper geometric database of the Vandermonde, GauB, Kummer, Saalschiitz, Dixon, 
Clausen and Dougall identities (see [^), and many more identities, can be obtained using 
sumrecursion. As examples, we consider the difficult cases of Clausen and Dougall:^ 

^Note that with REDUCE Version 3.5 we use the global operator summ instead of sum to denote the sum. 
*Note that the latter may need a large amount of computing time. 
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13 : summand : =f actorial (a+k-1) *f actorial (b+k-1) / 
(factorial (k) *f actorial (-1/2+a+b+k) ) * 
factorial (a+n-k-1) *f actorial (b+n-k-1) / 
(factorial (n-k) *f actorial (-1/2+a+b+n-k) ) $ 

14: suiiirecursion(summand,k,n) ; 

(2*a + 2*b + 2*n - l)*(2*a + 2*b + n - l)*sum(n)*n 

- 2*(2*a + n - l)*(a + b + n - l)*(2*b + n - l)*suin(n - 1) 

15 : summand : =pochhaiiimer (d , k) *poclihammer ( 1+d/ 2 , k) * 
pochhammer (d+b-a , k) *pochhammer (d+c-a , k) * 
pochhammer ( 1+a-b-c , k) *pochhammer (n+a, k) * 
pochhammer (-n , k) / (factorial (k) *pochhammer (d/2 , k) * 
pochhammer ( 1+a-b , k) *pochhammer ( 1+a-c , k) * 

pochhammer (b+c+d-a , k) *pochhammer ( l+d-a-n , k) *pochhammer ( l+d+n , k) ) $ 
16: sumrecursion (summand, k, n) ; 

(2*a - b - c - d + n)*(b + n - l)*(c + n - l)*(d + n)*sum(n - 1) 
+ (a-b-c-d-n+ l)*(a - b + n)*(a - c + n)*(a - d + n - l)*sum(n) 

corresponding to the statements 

\ {2a)nia + b)ni2b)n 



(a ,b ,1/2 — a — b — n , —n 
l/2 + a + b,l-a-n,l-b-n 



and 



(2a + 2b)n (a)n (b), 

^ . d,l + d/2,d+b — a,d + c — a,l + a — b — c,n + a, —n 
^ ^\d/2 ,1 -\- a — b ,1 -\- a — c ,b -\- c-\- d — a ,1 -\- d — a — n ,1 -\- d -\- n 

{d + 1)„ {b)n {C)n {l + 2a-b-C-d)r 



(a - d)nil + a-b)n{l + a-c)n{b + c + d- a)n 
(compare next section), respectively. 

Other apphcations of the Zeilberger algorithm are connected with the verification of identities. 
To prove the identity 

e. g., we may prove that both sums satisfy the same recurrence equation 

17: sumrecursion (binomial (n,k) "3, k,n) ; 

2 2 2 

8*(n - 1) *sum(n - 2) - sum(n)*n + (7*n - 7*n + 2)*sum(n - 1) 



6 



18: sumrecursion (binomial (n,k) '■2*binoiiiial(2*k,n) ,k,n) ; 



2 



2 



2 



8*(n - 1) *sum(n - 2) - suni(n)*n + (7*n - 7*n + 2)*sum(ii - 1) 

and finally check the initial conditions: 

19: sum(sub(ii=0, binomial (ii,k)"3) ,k,0,0) ; 

1 

20: sum(sub(n=0, binomial (n,k)"2*binomial(2*k,n)) ,k,0,0) ; 
1 

21 : sum(sub(n=l, binomial (n,k)"3) ,k,0, 1) ; 
2 

22: sum(sub(n=l, binomial (n,k)"2*binomial(2*k,n)) ,k,0, 1) ; 
2 



5 REDUCE operator HYPERRECURSION 



Sums for which the Zeilberger algorithm applies, in general are special cases of the generalized 
hypergeometric function 



with upper parameters {ai, 02, ... , ap}, and lower parameters {61, 62, ... , bg}. If a recursion 
for a generalized hypergeometric function is to be established, you can use the following 
REDUCE operator: 

• hyperrecursionCupper, lower ,x,n) determines a holonomic recurrence equation with 



upper parameters, and lower= {61, 62, • • • , bg} is the list of lower parameters depending 
on n. 

• hyperrecursionCupper , lower, x,n,j ) (j G IN) searches only for a holonomic recur- 
rence equation of order j. 





Therefore 
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23 : hyperrecursion({-n,b} , {c}, 1 ,n) ; 



(n - 1 + c - b)*sum(n - 1) - (n - 1 + c)*sum(n) 



establishes the Vandermonde identity 



2-^1 



—n 



c 



b 



1 



) 



(c-b) 

(C)n 



n 



whereas 



24 : hyperrecursion({d, l+d/2 ,d+b-a,d+c-a, 1+a-b-c ,n+a,-n} , 

{d/2 , 1+a-b , 1+a-c , b+c+d-a , 1+d-a-n , 1+d+n} , 1 , n) ; 



(2*a - b - c - d + n)*(b + n - l)*(c + n - l)*(d + n)*sum(n - 1) 
+ (a-b-c-d-n + l)*(a - b + n)*(a - c + n)*(a - d + n - l)*sum(ii) 

proves Dougall's identity, again. 

If a hypergeometric expression is given in hypergeometric notation, then the use of 
hyperrecursion is more natural than the use of sumrecursion. 
Moreover you may use the REDUCE operator 

• hyperterm (upper, lower, x,k) that yields the hypergeometric term 



with upper parameters upper= {ai, 02, . . . , flp}, and lower parameters lower= 
{61,62, • • • ,6J 

in connection with hypergeometric terms. 

The operator sumrecursion can also be used to obtain three-term recurrence equations for 
systems of orthogonal polynomials with the aid of known hypergeometric representations. 
By ([^, (2.7.11a)), the discrete Krawtchouk polynomials k^\x,N) have the hypergeometric 
representation 



and therefore we declare 

25: krawtchoukterin:=(-l)"n*p"n*binomial(NN,n)*hyperterm({-n,-x},{-NN}, l/p,k)$ 

and get the three three-term recurrence equations 



(Qi)fc • (Q2)fc • • • (Qp)fc 



{bi)k ■ {b2)k ■ ■ ■ {bq)k k\ 
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26: sumrecursioii(krawtchoukterni,k,ii) ; 

(x + 1 - 2*p - im*p + (2*p - l)*n)*sum(n - 1) 

- ((n - nn - 2)*(p - l)*sum(ii - 2)*p + suin(ii)*ii) 
27: suinrecursion(krawtchoukterm,k,x) ; 

- ((x - 1 + nn*p - n - 2*(x - l)*p)*suni(x - 1) + (x - 1 - nn) *siani(x) *p 

+ (p - l)*(x - l)*suin(x - 2)) 
28: sumrecursion(krawtchoukterni,k,NN) ; 
(x + 1 + n + (p - 2)*nn)*sum(nn - 1) 

- ((x + 1 - nn)*suin(nn - 2) - (n - nn)*(p - l)*suni(nn)) 
with respect to the parameters n, x, and N respectively. 

6 Simplification Operators 

For the decision that an expression at is a hypergeomctric term, it is necessary to find out 
whether or not a^/ak-i is a rational function with respect to k. For the purpose to decide 
whether or not an expression involving powers, factorials, F function terms, binomial co- 
efficients, and Pochhammer symbols is a hypergeomctric term, the following simplification 
operators can be used: 

• simplif y_gaimna(f ) simplifies an expression f involving only rational, powers and 
F function terms according to a recursive application of the simplification rule F(a+1) = 
a F (a) to the expression tree. Since all F arguments with integer difference are trans- 
formed, this gives a decision procedure for rationality for integer-linear F term product 
ratios. 

• simplify_coinbinatorial(f ) simplifies an expression f involving powers, factorials, 

F function terms, binomial coefficients, and Pochhammer symbols by converting factori- 
als, binomial coefficients, and Pochhammer symbols into F function terms, and applying 
simplif y.gamma to its result. If the output is not rational, it is given in terms of F 
functions. If you prefer factorials you may use 

• gammatof actor ial (rule) converting F function terms in factorials using T{x) — > (x— 1)!. 

The use of siniplify_combinatorial(f ) is a safe way to decide the rationality for any ratio 
of products of powers, factorials, F function terms, binomial coefficients, and Pochhammer 
symbols. 
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Example: 



29 : simplif y_combiiiatorial (sub (k=k+l , krawtchoukterm) /krawtchoukterm) ; 
(k - n)*(k - x) 



(k - nn)*(k + l)*p 

From this calculation, we see again that the upper parameters of the hypergeometric repre- 
sentation of the Krawtchouk polynomials are given by {— n, —x}, its lower parameter is {—n}, 
and the argument of the hypergeometric function is 1/p. 
Another example is 

30: simplif y_conibinatorial (binomial (n,k)) ; 



gamma(k + l)*gamma(n + 1 - k) 
31: ws where gammatof actor ial; 
factorial (n) 
factorial ( - (k - n) ) *f actorial (k) 

7 Tracing 

If you set 

32: on zb_trace; 

tracing is enabled, and you get intermediate results, see 0. 
Example for the Gosper algorithm: 

33: gosper (pochhammer (k-n,n) ,k) ; 



gamma (n + 1) 



k - 



1 



a(k)/a(k-l) := 



k - 



n 



- 1 



Gosper algorithm applicable 



p:= 1 



q 



:= k - 



1 
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r:= - (n + 1 - k) 
degreeboTind := 
1 

f:= 

n + 1 

Gosper algorithm successful 
pochhammer (k - n,n)*k 



n + 1 

Example for the Zeilberger algorithm: 
34: sumrecursion (binomial (n,k) '■2,k,n) ; 

2 

n 

F(n,k)/F(n-l,k) := 

2 2 
k - 2*k*n + n 

2 2 
k - 2*k*n - 2*k + n + 2*n + 1 

F(n,k)/F(n,k-1):= 

2 

k 

Zeilberger algorithm applicable 

applying Zeilberger algorithm for order := 1 

2 2 2 

p:= zb_sigma(l)*k - 2*zb_sigma(l)*k*n + zb_sigma(l)*n + n 

2 2 
q:= k - 2*k*n - 2*k + n + 2*n + 1 

2 

r:= k 

degreebound := 1 

2*k - 3*n + 2 
f .= 

n 
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2 2 2 3 2 

- 4*k *n + 2*k + 8*k*n - 4*k*n - 3*n + 2*n 
p.= 

n 

Zeilberger algorithm successful 

4*suni(n - l)*n - 2*suiii(n - 1) - suiii(n)*n 

35: off zb_trace; 

8 Global Variables and Switches 

The following global variables and switches can be used in connection with the ZEILBERG 
package: 

• zb_trace, switch; default setting on. Turns tracing on and off. 

• zb_direction, variable; settings: down, up; default setting down. 

In the case of the Gosper algorithm, either a downward or a forward antidifference is 
calculated, i. e., gosper finds gk with either 

O'k — gk — gk-l or Qk — Qk+l — Qki 

respectively. 

In the case of the Zeilberger algorithm, either a downward or an upward recurrence 
equation is given. Example: 

36: zb_direction:=up$ 

37: sumrecursion(binomial(ii,k) "2,k,ii) ; 

sum(n + l)*n + sumdi + 1) - 4*sum(n)*n - 2*sum(n) 

38: zb_directioii:=down$ 

• zb_order, variable; settings: any nonnegative integer; default setting 5. Gives the 
maximal order for the recurrence equation that sumrecursion searches for. 

• zb_factor, switch; default setting on. If off, the factorization of the output usually 
producing nicer results is suppressed. 

• zb_proof , switch; default setting off. If on, then several intermediate results are stored 
in global variables: 
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gosper_representation, variable; default setting nil. 

If a gosper command is issued, and if the Gosper algorithm is applicable, then the 
variable gosper_representation is set to the list of polynomials (with respect to k) 
{p,q,r,f } corresponding to the representation 

Ofc Pk <ik Qk+l J. 

Ok — Jk flfc , 



see [jl|. Examples: 

39: on zb_proof ; 

40: gosper(k*factorial(k) ,k) ; 

(k + l)*factorial(k) 

41: gosper_represeiitation; 

{k,k,l,l} 

42: gosper(l/(k+l)*biiiomial(2*k,k)/(n-k+l)*binoniial(2*n-2*k,n-k) ,k) ; 
(2*k - n + l)*(2*k + l)*binomial( - 2*(k - n) , - (k - n))*binoniial(2*k,k) 

(k + l)*(n + 2)*(n + 1) 

43: gosper_representation; 
{1. 

(2*k - l)*(k - n - 2) , 
(2*k - 2*n - l)*(k + 1) , 
n - 1 - 2*k 



} 

(n + 2)*(n + 1) 



zeilberger_representation, variable; default setting nil. 

If a sumrecursion command is issued, and if the Zeilberger algorithm is successful, then 
the variable zeilberger_representation is set to the final Gosper representation used, 
see PI. 
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9 Messages 



The following messages may occur: 

• ***** Gosper algorithm: no closed form solution exists 
Example input: 

gosper(f actorial(k) ,k). 

• ***** Gosper algorithm not applicable 

Example input: 

gosper(f actorial(k/2) ,k). 

The term ratio Uk/ak-i is not rational. 

• ***** illegal number of arguments 
Example input: 

gosper(k) . 

• ***** Zeilberger algorithm fails. Enlarge zb_order 
Example input: 

sumrecursion (binomial (n , k) *binomial (6*k , n) , k , n) 
For this example a setting zb_order:=6 is needed. 

• ***** Zeilberger algorithm not applicable 
Example input: 

sumrecursion (binomial (n/2,k) ,k,n) 

One of the term ratios f{n, k)/ f{n — 1, k) or f{n, k)/ f{n, k — 1) is not rational. 
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